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Abstract 


The  purpose  of  the  proposed  research  project  was  to  address  the  problem  of  com¬ 
munication  black-out  with  a  vehicle  during  a  sustained  hypersonic  flight.  Attempts  to 
communicate  through  dense  plasmas,  which  forms  around  an  aircraft  during  hyper¬ 
sonic  flight,  fail  when  the  transmitted  wave  frequency  is  below  the  plasma  frequency. 
Under  such  conditions,  the  transmitted  propagating  wave  becomes  evanescent  inside 
the  plasma  layer  and  its  amplitude  is  exponentially  attenuated.  It  was  demonstrated 
that  without  accounting  for  dissipation,  anomalous  plasma  transparency  to  radio¬ 
frequencies  can  be  achieved  via  amplification  of  evanescent  waves  using  a  diffraction 
grating,  periodic  plasma  modulation,  or  a  barrier.  For  dissipative  plasma,  adequate 
transmission  can  be  achieved  using  a  diffraction  grating  or  periodic  modulation.  Fur¬ 
thermore,  using  a  barrier,  reflectionless  transmission  and  absorption  can  be  achieved 
in  dissipative  plasma. 
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Status  of  Effort 


The  purpose  of  the  proposed  research  project  was  to  address  the  problem  of  com¬ 
munication  black-out  with  a  vehicle  during  a  sustained  hypersonic  flight.  Attempts  to 
communicate  through  dense  plasmas,  which  forms  around  an  aircraft  during  hyper¬ 
sonic  flight,  fail  when  the  transmitted  wave  frequency  is  below  the  plasma  frequency. 
Under  such  conditions,  the  transmitted  propagating  wave  becomes  evanescent  inside 
the  plasma  layer  and  its  amplitude  is  exponentially  attenuated.  The  main  goal  was 
to  develop  a  theoretical  framework  for  achieving  anomalous  plasma  transparency  to 
radio-frequencies. 

The  main  idea  of  the  proposed  method  was  to  amplify  the  evanescent  wave  through 
resonance  with  the  surface  wave,  and  then  use  that  evanescent  wave  to  transmit  a 
GPS  signal  through  dense  plasma.  In  order  to  achieve  resonance,  it  was  proposed  to 
use  periodic  plasma  modulation,  induced  by  a  diffraction  grating.  Similar  methods 
have  recently  been  used  in  near  field  optics  and  nano-photonics.  It  seemed  natural 
to  ask  whether  such  an  approach  can  be  used  to  solve  the  communications  problem. 

During  the  first  period  of  the  research  project,  we  concentrated  on  formulating 
mathematical  models  for  the  communication  problem  which  take  relevant  plasma 
properties  under  consideration.  We  worked  in  three  different,  but  related  directions: 

•  formulating  and  studying  the  transmission  problem  for  homogeneous  plasmas; 

•  understanding  the  properties  and  geometry  of  the  non-homogeneous  bounded  plasma 
that  forms  during  a  hypersonic  flight; 

•  formulating  and  studying  the  transmission  problem  for  non-homogeneous  plasma. 

In  particular,  we  have  demonstrated  that  plasma  transparency  can  be  achieved 
through  periodic  plasma  modulation  for  plasma  parameters  which  occur  during  a 
hypersonic  flight. 

During  the  second  year  of  the  research  project,  we  investigated  various  resonance 
conditions  which  could  be  relevant  for  hypersonic  flights  and  further  our  understand¬ 
ing  of  signal  transmission  through  dense  plasma.  In  particular,  we  concentrated  on 
specific  problems  related  to  electromagnetic  wrave  transmission,  such  as 

•  resonant  transmission  in  multi-layer  dense  plasma  structures; 

•  standing  wave  resonance  in  a  multi-layer  structure; 

•  resonant  transmission  through  dense  plasma  and  a  partially  modulated  boundary 
layer; 
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•  resonant  modes  and  tunneling  in  multi-layer  structures; 

•  effects  of  electron-atom  collisions; 

•  modeling  two-dimensional  bounded  plasmas  for  intermediate  and  high  pressures. 

In  particular,  we  have  made  the  first  attempt  to  achieve  resonant  transmission 
using  a  simplified  model  of  a  diffraction  grating. 

In  the  third  and  final  period  of  the  research  project,  we  have  concentrated  on 
signal  transmission  using  a  diffraction  grating.  In  particular,  we  formulated  a  general 
mathematical  model  which  employs  a  diffraction  grating.  During  this  period,  we  have 
been  working  in  the  following  directions: 

•  developing  a  model  which  describes  signal  transmission  through  dense  plasma  using 
a  diffraction  grating; 

•  understanding  the  parameters  of  the  diffraction  grating  which  are  needed  to  achieve 
plasma  transparency; 

•  understanding  plasma  transparency  for  different  incidence  angles; 

•  understanding  resonance  induced  by  periodic  plasma  modulation; 

•  understanding  resonance  induced  by  a  barrier; 

•  effects  of  dissipation  on  signal  transmission; 

•  modeling  two-dimensional  bounded  collisional  plasma; 

•  understanding  effects  of  oblique  magnetic  fields. 

The  research  has  been  done  in  collaboration  with  Andrei  Smolyakov,  Department 
of  Physics,  University  of  Saskatchewan,  Canada  and  Valery  Godyak,  RF  Plasma 
Consulting,  Brookline,  MA. 


Publications 

1.  N.  Sternberg  and  V.  Godyak,  The  Bohrn  plasma-sheath  model  and  the  Bohm 
criterion  revisited,  Special  Issue,  IEEE  Trans.  Plasma  Sci.,  35,  5,  1341,  2007. 

2.  V.  Godyak  and  N.  Sternberg,  Two-dimensional  cylindrical  plasma  at  low  gas 
pressure,  Plasma  Sources  Sci.  and  Techn.,  17,  025004,  2008. 

3.  N.  Sternberg,  C.  Sataline  and  V.  Godyak,  Influence  of  Ramsauer  effect  on  bounded 
plasmas  in  magnetic  fields,  35th  EPS  Conference  on  Plasma  Phys.  Hersonissos,  9-13 
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June  19,  EGA  Vol.  32,  P.-1.183,  2008. 


4.  N.  Sternberg  and  A.  I.  Smolyakov,  Resonant  transmission  of  electromagnetic  waves 
in  multilayer  dense  plasma  structures,  IEEE  Trans.  Plasma  Sci.,  37,  7,  2009. 

5.  E.  Fourkal,  I.  Velchev  and  A.  Smolyakov,  Energy  and  information  flow  in  super- 
lensing,  Plivs.  Rev.  A  79,  033846  2009. 

6.  A.  M.  Froese,  A.  I.  Smolyakov  and  D.  Sydorenko,  Nonlinear  skin  effect  in  a 
collisionless  plasma,  Phys.  of  Plasmas  16,  080704  2009. 

7.  N.  Sternberg  and  A.  Smolyakov,  Resonant  transmission  through  dense  plasmas 
via  amplification  of  evanescent  mode,  PIERS  Proceedings,  Moscow,  Russia,  August 
18-21,  2009,  p.1618. 

8.  N.  Sternberg  and  A.  I.  Smolyakov,  Resonant  transmission  through  dense  plasmas 
via  amplification  of  evanescent  mode,  PIERS  ONLINE,  vol.  5,  no.  8,  2009,  p.  781. 

9.  N.  Sternberg  and  A.  Smolyakov,  Resonant  transparency  of  a  three-layer  structure 
containing  the  dense  plasma  region,  Progress  In  Electromagnetics  Research,  PIER 
99,  37-52,  2009. 

10.  A.  Smolyakov,  S.  I.  Krashennikov,  E.  Furkal  and  N.  Sternberg,  Resonant  modes 
and  tunneling  in  multi-layer  metal  dielectric  structures  (to  be  submitted). 

Work  in  Progress 

1.  N.  Sternberg  and  A.  Smolyakov,  Resonant  transmission  of  electromagnetic  waves 
using  a  diffraction  grating. 

2.  A.  Smolyakov  and  N.  Sternberg,  Resonance  transmission  of  electromagnetic  waves 
induced  by  plasma  modulation. 

3.  N.  Sternberg  and  A.  Smolyakov,  Reflectionless  transmission  and  absorption  in 
dense  plasma  with  a  barrier. 

4.  N.  Sternberg  and  A.  Smolyakov,  Effects  of  dissipation  of  electromagnetic  wave 
transmission  through  dense  plasma. 

5.  N.  Sternberg  and  V.  Godyak,  Two-dimensional  collisional  discharge  plasma. 

6.  N.  Sternberg  and  J.-L.  Raimbault,  Effects  of  oblique  magnetic  fields  on  plasmas. 
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Conference  Presentations 


1.  “Anomalous  transparency  of  opaque  materials”,  Canadian  Association  of  Physi¬ 
cists  Congress,  June  18-20,  2007  Saskatoon,  Canada, 

2.  “The  Bohm  plasma-sheath  model  and  the  Bohm  criterion  revisited” ,  60  Gaseous 
Electronics  Conference,  Washington,  DC,  October  2-5,  2007. 

3.  “Two-dimensional  collisionless  weakly  ionized  plasma  in  fluid  approximation”,  60 
Gaseous  Electronics  Conference,  Washington,  DC,  October  2-5,  2007 

4.  “Influence  of  Ramsauer  effect  on  bounded  plasmas  in  magnetic  fields” ,  35  European 
Plasma  Conference,  June  13-June  19,  2008,  Crete,  Greece. 

5.  “Two-dimensional  gas  discharge  plasma”,  2009  Workshop  on  Radio-Frequency 
Discharges,  La  Presqu'il  de  Giens,  France,  17-20  May,  2009. 

6.  “Resonant  transparency  of  multi-layer  structures  with  opaque  materials”,  9th 
International  Conference  on  Mathematical  and  Numerical  Aspects  of  Wave  Propaga¬ 
tion,  15-19  June,  2009,  Pau.  France. 

7.  “Resonant  transmission  through  dense  plasmas  via  amplification  of  evanescent 
mode”,  Progress  in  Electromagnetics  Research  Symposium  (PIERS)  2009,  Moscow, 
Russia,  18-21  August,  2009. 

Presentations  at  the  AFOSR  T&E  Portfolio  Review 

1.  “Communication  through  plasma”,  Washington,  DC,  September  18-19,  2007. 

2.  “Communication  through  plasma  during  hypersonic  flight”,  Reston,  VI,  August 
5-6,  2008. 

3.  “Communication  through  plasma  during  hypersonic  flight”,  Dayton,  OH,  Septem¬ 
ber  21-23,  2009. 
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Main  Results 


1 .  Some  theoretical  aspects  of  resonant  transmission  induced  by  a 
diffraction  grating . 


Resonant  transmission  induced  by  a  diffraction  grating  was  first  studied  in  [1], 
[2].  There,  dense  plasma  was  bounded  by  a  diffraction  grating  on  each  side.  Such 
symmetric  configuration,  however,  is  not  applicable  for  conditions  that  occur  during 
hypersonic  flight,  where  one  could  expect  to  use  only  one  diffraction  grating  at  the 
antenna. 

In  our  model,  we  considered  the  configuration  shown  in  Fig.  1.  A  dense  plasma 
layer  of  width  l  =  b  — a  is  nested  between  a  semi-infinite  vacuum  layer  and  a  dielectric 
layer  of  positive  permittivity.  A  diffraction  grating  is  placed  inside  the  dielectric  layer 
at  a  distance  a  from  the  plasma  layer.  We  consider  the  regions  to  be  in  the  (y,  z)-plane, 
and  we  set  z  =  0  at  the  position  of  the  diffraction  grating.  In  this  way,  we  obtain 
four  distinct  regions:  the  semi-infinite  vacuum  region  V  to  the  left  of  the  plasma  layer 
(z  <  —6),  the  plasma  layer  P  (—6  <  z  <  —a),  the  dielectric  region  D\  (—a  <  z  <  0) 
to  the  left  of  the  diffraction  grating,  and  the  dielectric  layer  D 2  (z  >  0)  of  an  arbitrary 
width  to  the  right  of  the  diffraction  grating.  For  simplicity  of  exposition,  we  assume 
that  the  dielectric  regions  D\  and  Do  have  the  same  permittivity. 

In  a  general  medium,  the  dispersion  relation  of  the  incidence  wave  of  frequency  u 


H  =  Ho  exp (ikzz  +  ikyy  —  iut ) 


(i) 


is  given  by 


(2) 


where  c  is  the  speed  of  light.  For  vacuum,  6=1,  for  dense  plasma  e  =  ep  <  0,  and  for 
the  dielectric  material  e  =  €d  >  1.  Since  ky  =  (<j/c)  sin(0)  <  u/c,  as  long  as  e  >  1,  it 
holds  that  kz  >  0,  and  the  wave  is  propagating.  In  the  plasma,  however,  k2z  <  0,  the 
wave  is  evanescent  and  decays  rapidly.  This  is  illustrated  in  Fig.  2. 

Although  the  wave  in  the  plasma  region  is  decaying,  two  modes  exist  in  that  region, 
one  exponentially  decaying  and  the  other  one  exponentially  growing.  It  was  shown 
in  [3]  that  if  both  exponential  modes  are  combined  they  carry  energy  and  can  be 
used  for  signal  transmission.  The  question  becomes  how  to  increase  the  amplitude 
of  the  transmitted  incident  wave  inside  the  dense  plasma  layer.  It  is  well  known 
that  this  can  be  achieved  using  surface  wave  resonance.  Surface  waves  form  at  a 
sharp  interfaces  between  a  medium  with  positive  permittivity  and  a  medium  with 
negative  permittivity.  The  dispersion  relation  of  the  surface  wave  between  two  media 
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of  permittivity  ?i  and  e2  is  given  by 


7(fi)  _  7(^2) 

ei  e2 


where 


which  yields  the  dispersion  relation 

2  cj2  eie2 

y  c2  ei  +  e2 


(3) 

(4) 

(5) 


If  a  wave  is  in  resonance  with  the  surface  wrave,  then  the  wrave  vector  ky  for  that  wave 
has  to  satisfy  (5).  Obviously,  it  is  not  possible  to  achieve  surface  wrave  resonance  at 
the  plasma-vacuum  or  plasma-dielectric  interfaces,  and  the  incident  signal  cannot  be 
amplified  as  it  passes  through  those  different  regions. 

By  placing  a  diffraction  grating  inside  the  dielectric  region,  two  sidebands  are 
generated: 

H+  =  H+(z)  exp (i(ky  +  q)y  -  iuot)  (6) 


and 


H  =  H-(z)  exp (i(ky  —  q)y  —  iut)) 


(7) 


If  at  least  one  of  the  sidebands  is  in  resonance  with  the  surface  wave,  one  of  the 
following  dispersion  relations  is  satisfied: 


(ky  +  q?  =  % 


or 


(ky  -  qf  = 


U) 


ep  +  ed 
^p^d 


c2  tv  +  fd 


(8) 


(9) 


This  yields  the  value  of  the  diffraction  grating  wave  vector,  q ,  which  is  needed  to 
achieve  resonance  between  the  sidebands  and  the  surface  wave. 


The  sidebands  and  the  incidence  wave  are  independent.  Each  satisfies  the  wave 
equation.  Passing  through  the  region,  each  wrave  satisfies  continuity  conditions  at  the 
corresponding  interface.  At  the  diffraction  grating,  however,  all  waves  are  coupled 
through  the  corresponding  boundary  condition.  If  the  amplitude  of  the  sidebands 
becomes  large  due  to  surface  wave  resonance,  then  the  amplitude  of  the  incident 
wave  will  be  amplified.  Thus,  the  incidence  wave  is  not  amplified  directly  by  the 
surface  wave.  The  surface  wave  amplifies  the  sidebands  through  resonance,  and  the 
sidebands  amplify  the  incidence  wave. 
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We  have  found  that  the  theory  developed  above  can  be  applied  to  plasmas  under 
hypersonic  flight  conditions.  This  is  illustrated  in  Fig.  3  when  the  plasma  dissipation 
is  neglected.  For  the  plasma  parameters,  we  have  chosen  the  values  suggested  in  [4] 
and  [5]:  the  plasma  width  is  l  =  0.02m  and  the  plasma  frequency  is  ujpe/ (2n)  =  6GHz. 
We  chose  the  incidence  frequency  within  the  GPS  range:  o;/(27t)  =  1GHz.  Then, 
neglecting  dissipation,  the  plasma  permittivity  is  ep  =  —35. 

Our  preliminary  results  were  presented  at  two  conferences: 

(i)  9th  International  Conference  on  Mathematical  and  Numerical  Aspects  of  Wave 
Propagation,  15-19  June,  2009,  Pau,  France. 

(ii)  Progress  in  Electromagnetics  Research  Symposium  (PIERS)  2009,  Moscow, 
Russia,  18-21  August,  2009. 

Some  of  our  preliminary  results  can  be  found  in  the  following  publications: 

(i)  N.  Sternberg  and  A.  Smolyakov,  Resonant  transmission  through  dense  plasmas 
via  amplification  of  evanescent  mode,  PIERS  Proceedings,  Moscow,  Russia,  August 
18-21,  2009,  p.1618. 

(ii)  N.  Sternberg  and  A.  I.  Smolyakov,  Resonant  transmission  through  dense  plas¬ 
mas  via  amplification  of  evanescent  mode,  PIERS  ONLINE,  vol.  5,  no.  8,  2009, 
p.  781. 

2.  Modeling  wave  transmission  through  dense  plasma  with  a  diffrac¬ 
tion  grating . 

As  in  the  previous  section,  we  consider  the  configuration  shown  in  Fig.  1.  Sup¬ 
pose  an  electromagnetic  incident  wave  propagates  in  the  vacuum  region  V  toward 
the  plasma  layer  with  the  incidence  angle  6.  The  question  is  whether  in  the  given 
configuration,  a  modulation  by  the  diffraction  grating  can  result  in  resonant  trans¬ 
mission  of  the  incident  wave  through  the  dense  plasma  layer  if  the  plasma  frequency 
c jjpe  is  above  the  wave  frequency  cj.  To  answer  this  question,  we  will  use  the  Fourier 
approximation  only  up  to  the  first  harmonics. 

We  consider  the  p-polarization:  II  =  (//x,  0,0),  and  E  =  (0,  Ey,Ez).  In  each 
region,  the  incident  wave  is  periodic  in  the  y-direction,  but  becomes  evanescent  in 
the  ^-direction  when  it  enters  the  dense  plasma  layer.  Thus,  we  can  represent  the 
incident  magnetic  wave  in  the  form  Ho(z)elkyV~lu}t  with  ky  =  (w/c)  sin(0).  As  a 
result  of  the  modulation  by  the  diffraction  grating,  two  sidebands  are  induced,  both 
periodic  in  the  y-direction,  one  with  the  wave  vector  ky  4-  y,  the  other  with  the  wave 
vector  ky  —  y,  where  y  is  the  wave  vector  of  the  diffraction  grating.  In  order  to 
achieve  resonant  amplification  of  the  amplitude  of  the  evanescent  incident  wave,  the 
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sidebands  have  to  be  evanescent  in  the  2-direction.  We  therefore  represent  them  in 
the  form  H+(z)e^ky+q^v~lu;t  and  H_{z)el^ky~q^y~lu)t .  Altogether,  we  are  looking  for  the 
magnetic  wave  field  of  the  form 

Hx(y,  z,  t )  =  [Ho(z)eikyV  +  H+(z)ei{ky+q)y  +  H-(z)ei{ky-q)y]e~iult  (10) 


as  it  is  transmitted  through  all  regions. 

Recall  that  for  a  double  non-homogenous  medium  (in  y  and  in  z)  the  Maxwell 
equations  yield  the  following  wave  equation: 


d  (l dHA  d  (ldHx\ 
£ dy  \£  dy  j  £ dz  \e  dz  ) 


U)  rT 

H — =  0 
c1 


(11) 


where  e  is  the  permittivity.  We  assume  that  e  is  defined  by  the  following  step- function: 
e(y,  z)  =  1  if  (j/,2)  is  in  the  vacuum  region,  e(y,  z)  =  Ed  if  (y,  z)  is  in  the  dielectric 
region,  and  e(y,  z)  =  ep  —  1  —  ^/(c o2  +  iveuo)  if  (y,  z)  is  in  the  plasma  region.  Here, 
ve  is  the  electron-atom  collision  frequency.  Thus,  each  of  the  regions  is  homogeneous, 
and  there  is  a  sharp  interface  between  them.  Note  that  ep  <  0  for  a  dense  plasma 
where  ujpe  >  uo  and  ve  =  0. 

In  order  to  account  for  the  diffraction  grating,  we  modify  (11)  as  follows: 

d  (ldHx\  d  (ldHx\  u2  TT  J2  7  ,  ,  NNC/,rr 

%  \7 ~ai )  +  \7 17 j  +  +  ^  (£» +  “cos  {qv))  1  (2)  H‘  =  0  (12) 


where  a  is  the  modulation  parameter,  eg  is  the  permittivity  of  the  diffraction  grating 
and  hg  is  its  width,  5  is  the  5-function.  Note  that  outside  the  diffraction  grating, 
equations  (11)  arid  (12)  coincide.  Equation  (12)  is  similar  to  the  one  found  in  [1],  [2], 
but  it  is  somewhat  more  general. 

At  the  plasma-vacuum  interface  and  at  the  plasma-dielectric  interface,  the  Maxwell 
equations  yield  continuity  of  Hx  and  of  e~ldHx/ dz,  i.e.,  at  z  —  —b  and  z  =  —a, 

[Hx]t  =  0  and  [1ft  =  °  (13) 

where  “+”  and  ”  indicate  the  corresponding  limits  in  the  z  direction  from  the 
right  and  from  the  left  of  the  interface  for  all  y.  At  the  diffraction  grating,  Hx  is 
continuous,  and  integrating  equation  (12)  yields 


where 
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In  order  to  find  the  corresponding  wave  equations  for  the  incident  wave  and  the 
sidebands  in  each  region,  we  substitute  (10)  into  (11)  and  obtain: 


^  -  k2yH0  +  ^eH0  =  0 

(16) 

dz2  (ki y  +  9)  H+  +  c2  fH+  -  0 

(17) 

( PH_  „  x2„  oj2  TT 

dz ?  <fc»  «>  H-  +  c‘,H-=a 

(18) 

The  boundary  conditions  (13)  and  (14)  yield  the  following  boundary  conditions  for 
the  incident  w’ave  and  the  sidebands: 


[Ho]t  =  0  [H+]t  =  0  [ff_]t  =  0  (19) 


at  the  plasma-vacuum  interface,  the  plasma-dielectric  interface  and  at  the  diffraction 
grating.  Furthermore, 


,l<Mo]  + 

[ld/f+]+  = 

rl  dH_, 

U  dz  i- 

U  dz  i- 

U  dz  J 

at  the  plasma-vacuum  interface  and  the  plasma-dielectric  interface.  Moreover, 

[77]- =  ~  yh-(0)  <21) 

=  °)-Yw«(°)  (22> 

[^1- =  -k«H-m  ~  y^0)  <23) 

at  the  diffraction  grating.  Here,  as  above,  “+”  and  ”  indicate  the  corresponding 
limit  from  the  right  and  from  the  left  along  the  2  axis.  Note  that  the  incident  wave 
and  the  sidebands  are  coupled  through  the  boundary  conditions  at  the  diffraction 
grating. 

In  each  region  shown  in  Fig.  1,  the  general  solution  of  the  equations  (16),  (17) 
and  (18)  is  of  the  form 


Hq  —  >l(exp(7s)  +  T  exp(— 72)) 

(24) 

H+  =  >l+(exp(7+2)  +  F+exp(— 7+2)) 

(25) 

//_  =  J4_(exp(7_2)  +  r~  exp(— 7-2)) 

(26) 
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where 

72  =  kl  -  (27) 

(7+)2  =  (^y  +  Q)2  ~  (28) 

(7-)2  =  (ky  -  qf  -  (29) 

Note  that  A,  A+  and  A~  denote  the  amplitudes  of  the  corresponding  incident 
waves,  while  I\  F+  and  T-  are  the  corresponding  reflection  coefficients. 

The  behavior  of  Hq,  H+  and  //_  depends  on  each  region.  In  the  vacuum  region, 
£  =  1,  ky  <  uj/c,  and  the  wave  Hq  is  propagating.  Hence, 

7  =  hv  (30) 

where 

^  =  cos(0)  (31) 

In  the  plasma  region,  Ho  is  evanescent,  and 

7  =  — 7 p  (32) 

where  ( _ 

7P  =  =  ^\/sin2(0)  -ep  (33) 

The  transmitted  wave  is  propagating  in  the  dielectric  region.  Thus,  Ed  >  sin2(0), 
and 

7  =  »7rf  (34) 

w'here  _ 

7d  =  ^yed-sin^)  (35) 

As  mentioned  above,  in  order  to  induce  resonance  with  the  evanescent  incident 

wave,  the  sidebands  H+  and  H_  have  to  be  evanescent.  This  is  achieved  if  q  > 
ky  +  ( uj/c)y/Ed ■  Furthermore,  the  sidebands  are  decaying  away  from  the  diffraction 


grating.  Thus,  in  the  dielectric  region 

A, 

-> 

+ 

II 

and 

i 

ii 

s?l 

(36) 

in  the  dielectric  region  Di , 

+ 

II 

1 

£+ 

and 

7  =  ~ld 

(37) 
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in  the  vacuum  region  V, 


7+  =  lv 


and  in  the  plasma  region  P, 


7+  =  Ip 


where 


~fv  =  V  +  ?)2  _ 


7 d  =  V  +  9)2  - 


Ip  =  V  (fc2/  +  # 


ur 


X2 


and 

7 

=  % 

(38) 

and 

7" 

— 

(39) 

1, 

UJ2 

and 

%  = 

y  (kv  -  q)  ~ 

n 

(40) 

and 

7 d  = 

1 

CM 

1 

—£d 

cz 

(41) 

and 

7 p 

1 

CM 

1 

oo2 

2 

Cz 

(42) 

Note  that  for  a  normal  incident  wave,  0  =  0,  and  H+  =  H  ,  which  simplifies  the 
computations. 

General  conditions  for  resonant  signal  transmission  via  a  multy-layer  structure 
were  studied  in  [6],  [7],  [8].  In  those  papers,  the  transmission  problem  was  formulated 
by  the  method  of  transfer  functions,  where  each  uniform  region  is  represented  by  a 
2x2  matrix.  This  method  leads  to  a  cumbersome  and  computationally  intense  result, 
wThich  is  represented  in  the  form  of  a  product  of  all  of  these  matrices  and  their  inverses. 
Instead,  we  used  the  impedance  method  which  relates  the  refection  coefficient  to  the 
load  impedance  of  the  vacuum  region  located  to  the  left  of  the  plasma.  In  particular, 
if  that  load  impedance  equals  the  vacuum  impedance  then  the  refection  coefficient 
vanishes  and  total  transmission  occurs. 


The  impedance  method  has  been  used  in  the  literature  to  study  propagating  waves. 
We  have  adapted  this  method  to  include  evanescent  waves.  The  main  advantage  of 
the  impedance  method  is  in  its  simplicity.  It  uses  a  recursion  formula  to  find  the 
impedance  at  each  interface,  which  lends  itself  to  generalizations  to  any  number  of 
layers  and  to  non-homogeneous  plasmas.  We  have  first  applied  this  method  to  two- 
and  three-layer  structures.  It  enabled  us  to  obtain  analytical  and  numerical  results 
for  electromagnetic  wave  transmission  in  a  simple  and  transparent  way  and  gain  a 
better  understanding  of  that  process. 

Our  results  concerning  multi-layer  structures  and  a  detailed  discussion  of  the 
impedance  method  can  be  found  in  the  following  publications: 

(i)  N.  Sternberg  and  A.  I.  Smolyakov,  Resonant  transmission  of  electromagnetic 
waves  in  multilayer  dense  plasma  structures,  IEEE  Trans.  Plasma  Sci.,  37,  7,  2009. 
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(ii)  N.  Sternberg  and  A.  Smolyakov,  Resonant  transparency  of  a  three-layer  struc¬ 
ture  containing  the  dense  plasma  region,  Progress  In  Electromagnetics  Research, 
PIER  99,  37-52,  2009. 

In  an  attempt  to  understand  the  effects  of  the  diffraction  grating,  we  have  ap¬ 
plied  the  impedance  method  to  the  model  presented  above  and  obtained  analytical 
transparency  conditions  which  show  that  the  total  wave  transmission  can  indeed  be 
achieved  through  surface  wave  resonance  with  the  sidebands  for  an  appropriate  choice 
of  the  parameters  of  the  diffraction  grating.  Those  theoretical  results  are  being  pre¬ 
pared  for  publication. 

3.  Effects  of  the  incidence  angle  and  diffraction  grating  parameters  on 
resonant  wave  transmission  through  dense  plasma  without  dissipation 

Studying  the  diffraction  grating  model  described  in  the  previous  section,  we  have 
found  that  resonant  transparency  can  be  induced  if  the  appropriate  parameters  are 
controlled.  There  are  seven  parameters  of  the  problem:  the  dielectric  permittivity  q, 
the  diffraction  grating  permittivity  the  position  of  the  diffraction  grating  z  =  a, 
the  diffraction  grating  wave  vector  q ,  the  diffraction  grating  modulation  amplitude  a, 
the  diffraction  grating  width  hg  and  the  incidence  angle  6 .  Our  analysis  has  shown 
that  the  parameter  a  is  a  function  of  the  plasma  parameters,  the  incidence  angle  and 
the  dielectric  permittivity.  In  our  computations,  we  set  hg  =  0.05m,  and  we  assumed 
a  metal  diffraction  grating,  i.e.,  eg  >  0.  In  practice,  eg  is  fixed  and  a  <  eg.  For  a 
plasma  without  dissipation,  we  chose  eg  =  2  and  a  =  1.8. 

We  were  able  to  demonstrate  that  resonant  transparency  can  be  achieved  by  ad¬ 
justing  q  and  placing  the  diffraction  grating  appropriately.  This  is  illustrated  in  Fig.  4 
for  Cd  =  1.2  and  different  incidence  angles.  The  resonant  position  of  the  diffraction 
grating  depends  on  the  incidence  angle.  In  fact,  a  larger  incidence  angle  requires  a 
larger  resonant  value  of  a.  For  example,  a  —  0.06m  for  6  =  0,  a .  =  0.07m  for  6  =  7r/6, 
a  =  0.08m  for  6  —  n/ 4,  and  a  =  0.09m  for  0  =  7r/3. 

Fig.  5  shows  that  the  position  of  the  diffraction  grating  plays  a  significant  role 
in  resonant  signal  transmission.  Perturbations  in  the  parameter  a  can  significantly 
reduce  signal  transmission.  In  practice,  however,  the  exact  incidence  angle  is  not 
known,  and  therefore  the  effects  of  the  incidence  angle  need  to  be  reduced.  Fig. 
6  shows  that  this  can  be  done  by  increasing  the  dielectric  permittivity.  For  larger 
dielectric  permittivity  the  effect  of  the  distance  from  plasma  to  diffraction  grating  for 
different  incidence  angles  is  less  pronounced.  In  addition,  larger  values  of  require 
smaller  values  of  a  for  achieving  resonance,  which  may  be  important  in  application. 
Fig.  7  shows  that  adequate  signal  transmission  can  indeed  be  achieved  for  all  incidence 
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angles  by  choosing  an  average  value  of  a  and  a  sufficiently  large  e d- 

More  detailed  theoretical  results  are  being  prepared  for  publication. 

4 .  The  diffraction  grating  model  with  plasma  dissipation . 

All  our  results  so  far  neglected  the  electron-atom  collisions  in  the  plasma,  i.e., 
ve  =  0.  In  practice,  however,  electron-atom  collisions  cannot  be  neglected.  Figs. 

8  and  9  show  the  reflection,  transmission  and  absorption  coefficients  for  9  =  0  and 
0  =  7r/3  if  there  is  no  diffraction  grating  present.  As  can  be  seen,  for  small  values  of 
i/e,  most  of  the  signal  is  reflected  by  the  plasma,  some  is  absorbed,  and  very  little  is 
transmitted.  Indeed,  for  ve  =  0.01  and  9  =  0,  the  reflection  coefficient  is  R  =  99.4%, 
while  the  transmission  coefficient  is  about  T  =  0.3%,  and  the  absorption  coefficient  is 
about  A  =  0.3%.  This  behavior  is  similar  as  the  angle  is  increased.  For  example,  for 
ve  —  O.OIcj  and  9  =  7r/3,  we  have  found  R  =  98.4%,  T  =  0.8%,  and  A  =  0.8%.  For 
intermediate  values  of  ve  most  of  the  signal  is  absorbed  by  the  plasma.  Eventually, 
for  large  values  of  ue  much  of  the  signal  is  transmitted.  This  is  due  to  the  fact,  that 
when  plasma  dissipation  is  taken  into  account  ,  ep  becomes  complex,  and  as  shown  in 
Fig.  10  the  plasma  eventually  becomes  a  dielectric  for  sufficiently  large  values  of  i/e. 

Using  a  diffraction  grating  can  significantly  improve  signal  transmission,  although 
resonant  transmission  cannot  be  achieved.  For  our  computations  we  chose  td  =  1.2, 
eg  =  4,  a  =  3.8  and  a  =  0.062m.  Figs.  11-13  show  the  transmission,  reflection  and 
absorption  coefficients,  respectively,  for  ve  =  0.01a;.  In  this  case,  the  corresponding 
extreme  values  are  Tmax  =  43%,  Rmtn  =  14%  and  Amax  =  43%.  In  comparison,  for 

9  =  7r/3,  Tmax  =  50%,  Rmin  =  6%  and  Amax  =  44%.  Thus,  by  increasing  the  incidence 
angle,  the  reflection  coefficient  can  be  decreased  and  the  transmission  coefficient  can 
be  increased.  Observe  that  the  bandwidth  is  wider  when  dissipation  is  present,  but 
by  increasing  the  incidence  angle,  the  bandwidth  becomes  even  wider  (see  Fig.  14). 

It  is  interesting  to  note,  for  both,  9  =  0  and  9  =  7r/3  there  is  only  one  value  of  the 
diffraction  grating  wave  vector  q  which  will  produce  partial  signal  transmission.  We 
have  found  that  this  is  the  case  for  9  >  57r/21.  For  9  <  hn/21  two  such  values  of  q 
exist.  The  transition  from  one  resonant  value  of  q  to  two  such  values  takes  place  at 
about  9  =  57r/21.  As  the  incidence  angle  is  decreased,  the  two  values  of  q  are  moving 
closer  to  each  other,  and  eventually  coalesce  for  0  =  0.  This  is  illustrated  in  Figs. 
15-17. 

Our  results  on  dissipation  are  preliminary.  More  studies  are  needed  in  this  direc¬ 
tion. 


15 


5.  Transmission  of  electromagnetic  waves  through  dense  plasma  using 
a  barrier . 

In  order  to  model  transmission  of  electromagnetic  waves  using  a  barrier  instead  of 
a  diffraction  grating,  we  set  a  =  0  in  equation  (12).  Deriving  analytical  transparency 
condition,  we  were  able  to  show  that  using  this  simple  model  one  can  obtain  resonant 
transmission.  The  resonance  in  this  case  is  due  to  standing  waves  rather  than  to 
surface  modes.  We  have  studied  this  type  of  resonance  in  [9].  Accounting  for  plasma 
dissipation,  we  were  able  to  demonstrate  reflectionless  transmission  and  absorbtion 
of  the  incident  electromagnetic  wave. 

Our  results  are  being  prepared  for  publication. 

6.  Transmission  of  electromagnetic  waves  through  periodic  plasma 
modulation . 

The  idea  of  using  a  periodic  modulation  to  transmit  electromagnetic  waves  trough 
opaque  material  can  be  found  in  [10].  We  were  able  to  demonstrate  that  this  method 
can  be  adapted  to  transmit  a  signal  of  GPS  frequency  through  dense  plasma.  While 
in  [10],  the  whole  plasma  layer  was  modulated,  we  have  found  that  transparency  can 
also  be  achieved  by  partial  plasma  modulation. 

In  order  to  model  signal  transmission  through  dense  plasma  by  periodic  modula¬ 
tion,  we  simplified  equation  (12)  by  setting  eg  =  0.  We  are  currently  in  the  process  of 
developing  a  theoretical  framework  which  will  clarify  what  types  of  resonance  are  pos¬ 
sible  in  this  situation.  Our  goal  is  to  describe  the  relationship  between  the  modulation 
amplitude  and  the  modulation  wave  vector  which  are  needed  to  induce  resonance.  We 
intend  to  submitted  our  results  for  publication. 

7.  Modeling  two-dimensional  bounded  plasmas  for  intermediate  and 
high  pressures . 

In  order  to  successfully  address  the  communication  black-out  problem  during  hy¬ 
personic  fight,  one  needs  to  understand  the  plasma  that  forms  around  the  vehicle  in 
this  situation.  We  started  by  modeling  a  bounded  collsionless  cylindrical  plasma  and 
studied  the  behavior  of  the  plasma  parameters.  We  are  currently  modeling  a  bounded 
two-dimensional  cylindrical  plasma  which  takes  ion-atom  and  electron-atom  collisions 
into  account  [11].  Our  preliminary  results  were  presented  at  the  2009  Workshop  on 
Radio-Frequency  Discharges,  La  Presqu’il  de  Giens,  France,  17-20  May,  2009.  A 
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paper  on  two-dimensional  collisional  gas  discharge  plasma  is  being  prepared  for  pub¬ 
lication  in  collaboration  with  V.  Godyak. 

8.  Plasma  under  the  action  of  a  magnetic  field . 

It  is  possible  that  a  plasma  modulation  can  indeed  be  obtained  using  a  weak 
magnetic  field.  In  general,  it  is  suggested  in  the  literature  that  a  magnetic  field  can 
facilitate  signal  transmission.  For  practical  applications,  this  magnetic  field  has  to  be 
weak.  The  question  arises:  What  does  “weak”  really  mean? 

In  [12],  we  have  studied  a  cylindrical  plasma  with  an  axial  magnetic  field,  and 
found  that  sufficiently  weak  magnetic  fields  have  no  effect  on  plasmas  and  in  this  case 
the  Boltzmann  equilibrium  for  the  electrons  can  be  assumed.  For  our  computations, 
we  used  a  Ramsauer  gas.  In  [13],  we  extended  this  result  to  non-Ramsauer  gases  and 
found  simple  transformation  formulas  between  the  corresponding  results.  We  have 
found  that  in  both  cases,  for  Ramsauer  and  non-Ramsauer  gases,  if  the  plasma  is 
affected  by  a  strong  magnetic  field,  then  the  Boltzmann  equilibrium  is  violated.  We 
introduced  a  magnetic  parameter  that  gives  a  criterion  for  the  cut-off  between  the  two 
regions  of  weak  (practically  negligible)  and  strong  (dominant)  magnetic  field  effects. 

We  are  currently  conducting  similar  investigations  in  order  to  understand  the  ef¬ 
fects  of  an  oblique  magnetic  field  on  bounded  plasmas  and  the  transition  from  weak 
to  strong  magnetic  fields.  This  work  is  done  in  collaboration  with  J.-L.  Raimbault, 
Ecole  Polytechnique,  Paris,  France  and  will  be  published. 
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Figure  1:  Schematic  representation  of  the  multi-layer  problem  with  a  diffraction  grating. 
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Figure  2:  A  propagating  normal  incidence  wave  in  the  vacuum  layer  enters  the  dense  plasma  layer 
becoming  evanescent  (red)  and  continues  propagating  in  the  dielectric  layer  with  a  small  amplitude. 
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Figure  3:  Resonant  electromagnetic  wave  transmission  through  a  dense  plasma  layer  using  a  diffrac¬ 
tion  grating  for  normal  incidence. 
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Figure  4:  Transmission  coefficient  versus  diffraction  grating  wave  vector  for  =  1.2  and  for  different 
incidence  angles:  6  =  0  (blue),  6  =  7r/6  (red),  6  =  tt/4  (green),  and  6  =  7r/3  (purple). 
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Figure  5:  Transmission  coefficient  versus  diffraction  grating  wave  vector  for  ej  =  1.2,  a  =  0.06m,  and 
for  different  incidence  angles:  0  =  0  (blue),  0  =  7t/6  (red),  0  =  it /A  (green),  and  0  =  7t/3  (purple). 
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Figure  6:  Resonant  distance  from  the  plasma  layer  to  the  diffraction  grating  versus  dielectric  per¬ 
mittivity  for  different  incidence  angles:  6  =  0  (blue),  9  =  7r/6  (red),  6  =  7r/4  (green),  and  6  =  7r/3 
(purple). 


23 


1 


0.9 


0.8 


c  0.7 

.2 

I  0.6 

8 

1  0.5 
'«/) 

2 

E  0.4 

co 

c 


0.2 


0.1  - 

Ol  ^  1  '  '  - ^ - - 

80  85  90  95  100  105  110 

diffraction  grating  wave  vector 


Figure  7:  Transmission  coefficient  versus  diffraction  grating  wave  vector  for  e<*  =  10,  a  =  0.0172m, 
and  for  different  incidence  angles:  0  —  0  (blue),  0  =  7t/6  (red),  0  =  7t/4  (green),  and  6  =  7t/3 
(purple). 
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Figure  8:  Transmission  coefficient  (blue),  reflection  coefficient  (green)  and  absorption  coefficient 
(red)  versus  normalized  electron  collision  frequency  for  0  =  0  when  no  diffraction  grating  is  present. 
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Figure  9:  Transmission  coefficient  (blue),  reflection  coefficient  (green)  and  absorption  coefficient 
(red)  versus  normalized  electron  collision  frequency  for  6  =  7r/3  when  no  diffraction  grating  is 
present. 
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Figure  10:  Real  and  imaginary  part  of  the  plasma  permittivity  versus  the  normalized  electron 
collision  frequency. 
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Figure  11:  Transmission  coefficient  versus  the  diffraction  grating  wave  vector  for  ve  =  0.01a;  and 
0  =  0 . 
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Figure  12:  Reflection  coefficient  versus  the  diffraction  grating  wave  vector  for  ve  =  0.01  uj  and  0  —  0. 
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Figure  13:  Absorption  coefficient  versus  the  diffraction  grating  wave  vector  for  ut  —  0.01a;  and  0 
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Figure  14:  Transmission  coefficient  versus  the  diffraction  grating  wave  vector  for  ve  —  O.Olw  and 
9  =  7t/3. 
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Figure  15:  Transmission  coefficient  versus  the  diffraction  grating  wave  vector  for  ve  =  0.01a;  and 
0  =  5tt/21. 
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Figure  16:  Transmission  coefficient  versus  the  diffraction  grating  wave  vector  for  ve  —  0.01a;  and 
9  =  7t/6. 
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Figure  17:  Transmission  coefficient  versus  the  diffraction  grating  wave  vector  for  vt  —  0.01a;  and 
6  —  7t/48. 
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